New Chiral Universality Class in a Frustrated Three-Leg Spin Ladder 
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We study a model of three S = 1/2 antiferromagnetic Heisenberg spin chains weakly coupled by 
on-rung and plaquette-diagonal interchain interactions. It is shown that the model exhibits a critical 
phase with central charge C = 2 and belongs to the class of "chirally stabilized" liquids recently 
introduced by Andrei, Douglas, and Jerez. By allowing anisotropic interactions in spin space, we 
find an exact solution at a Toulouse point which captures all universal properties of the model, 
including the SU(2) symmetric case. At the new critical point the massless degrees of freedom are 
described in terms of an effective S =1/2 Heisenberg spin chain and two critical Ising models. We 
discuss the spectral properties of the model, compute spin-spin correlation functions and estimate 
the NMR relaxation rate. 
PACS No: 75.10.Jm, 75.40.Gb 



In close parallel with qualitative difference between 
the integer-spin and half-integer-spin antiferromagnetic 
chains, first predicted by Haldane the universal prop- 
erties of standard spin ladders, i.e. those with the inter- 
chain exchange interaction (Jj_) across the rungs, dra- 
matically depend on the parity of the number of legs 
0). Ladders with even number of legs are disordered 
spin liquids with a finite gap in the excitation spectrum, 
while in odd-legged ladders there exists a gapless branch 
in the spectrum implying that the spin correlations de- 
cay algebraically. In the low-energy limit, the critical 
behavior of the odd-legged ladders is characterized by 
the central charge C = 1 corresponding to an effective 
5 = 1/2 Heisenberg spin chain || . The interesting ques- 
tion is whether frustration can lead to new types of in- 
frared (IR) behavior not encountered by the above de- 
scribed two "even-odd" scenarios. The existence of real 
quasi- ID antiferromagnets with a zigzag interchain inter- 
action, such as Cs 2 CuCh [f§ and Cu2(C 5 H l2 N 2 )CU [§, 
indicates that this question is not purely academic. The 
role of frustration has recently been addressed in a two- 
chain model with a weak zigzag interchain coupling. For 
an isotropic antiferromagnetic interaction the spectrum 
has an exponentially small gap, and the ground state is 
spontaneously dimerized [|6||7|]. However, an easy plane 
XXZ-type anisotropy supports a gapless phase charac- 
terized by nonzero local spin currents polarized along the 
anisotropy axis and algebraically decaying incommensu- 
rate spin correlations fa]. 

Recently another type of frustrated two-leg ladders was 
discussed In addition to the standard on-rung cou- 
pling Jl, the new model also includes an interaction J x 
along both diagonals of elementary plaquettes. In this 
letter, we consider a three-chain generalization of such a 
model: 



Ti — J|| Sj,n ' Sjjn+l + J±So,n ■ (Sl,„ + S2,n) 

n j=0,l,2 

+ Jx [(Sl,ri + S2, n ) ' So,«+l + (Sl >n +l + S2, n +l) ' Sq,tJ} (1) 

where frustration introduced by the interaction J x shows 
up in a highly nontrivial way. We shall show that in a 
certain range of parameters the model displays a new 
type of critical behavior in the IR limit, with the central 
charge C — 2, rather than C — 1 as in usual (nonfrus- 
trated) three- leg ladders ]1C| ]. This critical behavior is 
identified as the universality class of "chirally stabilized" 
fluids, recently introduced by Andrei et al. [0. 

We shall assume that < Jj_, J x <C J\\. Under this 
condition a continuum description can be adopted in 
which the spin densities in each chain are represented as 
Sj (x) = Jj (x) + (-l) x l a ° rij (x), where Jj = Jjr + Jjl 
and Tdj are, respectively, the smooth and staggered parts 
of the magnetization. In the continuum limit the Hamil- 
tonian density of the original model (Q) takes the form: 



E • 

J'=0,l,2 



gi n ■ (ni + n 2 ) + g 2 Jo ■ (Ji + J2) 



(2) 



Here the first line describes three decoupled chains in 
terms of three critical SU[2)i Wess-Zumino-Novikov- 
Witten (WZNW) models (T^, Vj ~ Jya being the spin 
velocities. The second line in (Q) includes the inter- 
chain coupling terms with the constants: g\ = (Jj_ — 
2J x )ao, <? 2 = (Jj. + 2J x )ao- We stress that in the con- 
tinuum limit there is no marginally relevant twist per- 
turbation, iii ■ d x rij, which appears in the description of 
spin ladders with small zigzag interchain coupling ^,^3) . 

The two interaction terms in (^|) are of different na- 
ture. The first one is a relevant perturbation with scal- 
ing dimension d = 1. The second term describes an in- 
teraction between the total current of the surface chains, 
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I = Ji + J2 and that of the middle chain. This interac- 
tion is only marginal and, as long as g\ is not too small, 
can be discarded. As a result, for generic values of g\ and 
g2, the low energy physics of our model will be essentially 
that of the standard three-leg ladder, and frustration will 
plays no role (except for renormalization of the mass gaps 
and velocities [fL4|). The important point, though, is that 
in contrast with non-frustrated ladders, the two coupling 
constants 171 , gi can vary independently, and there exists 
a vicinity of the line J± — 2J X (g\ = 0) where the low- 
energy properties of the model are mainly determined by 
the current-current interchain interaction. Remarkably 
enough, exactly at g\ — an anisotropic version of this 
model is exactly solvable at a Toulouse point. The so- 
lution describes a new fixed point with a larger central 
charge C — 2 which, as we shall demonstrate, govern the 
physics for sufficiently small g\ . Our strategy will be first 
to present our solution and discuss the physical proper- 
ties of the system at the special point Jj_ = 2 J x , and 
then to explore its neighborhood for g\ small. 

Toulouse point solution. We start by neglecting 
the interaction between currents of the same chirality. 
Such terms lead to renormalization of the velocities of 
the excitations which can be effectively taken into ac- 
count by allowing the surface chain velocities [v\ = v%) 
to be different from that of the middle chain vq. In 
the noninteracting case (172 = 0), the central charge of 
two surface chains is C = 2. On the other hand, the 
total current I, being the sum of two SU(2)i currents, 
satisfies the SU(2j2 Kac-Moody algebra. Since the cen- 
tral charge of the SU(2)2-symmetric WZNW model is 
C = 3/2, some degrees of freedom should account for 
the missing C — 1/2. Those are associated with a dis- 
crete (Z2) nonmagnetic 1 «-> 2 interchange symmetry, 
and remain decoupled and critical with the central charge 
C = 1/2. The simplest way to see this is to exploit the 
equivalence SU(2) x SU(2) SO(4) and use the repre- 
sentation of two SU(2)i currents in terms of a quadru- 
plet of real (Majorana) fermions, £° and £ a (a = 1,2,3) 

0: IS = A (Jai-J a2 ) a = i where 

a = R, L. Then one easily finds that the contribution 
of the massless Majorana fermion £° decouples from the 
rest of the spectrum, and the effective Hamiltonian at 
the special point g\ = reads: 

n = -iy - e L d x e L )+H{i,u (3) 

All nontrivial physics is incorporated in the current 
dependent part of the Hamiltonian, Tt, describing 
marginally coupled SU(2) 2 and SU(2)i WZNW models. 
Notice that 7i in turn separates into two commuting and 
chirally asymmetric parts: TL = TL1+TL2, ([Hi,H2] = 0), 
where 

H\ = -y- Ir ■ Ir H ^ J L ■ Jol + 52 Ir • Jol, (4) 

with TL2 obtained from "Hi by inverting chiralities of all 
the currents. The Hamiltonian (^ resembles the two- 
channel Kondo model [jl5|: in the latter case the SU(2) 2 



current Ir describes two-flavor spin excitations of the 
conduction electrons, while the SU(2)i current Jol is 
replaced by the local spin density of the impurity spin 
S = 1/2. Moreover, one can show that (0) corresponds 
to the Hamiltonian of the two channel Kondo lattice away 
from half-filling with a nearest neighbour interaction be- 
tween the impurity spins. 

A simple renormalization group analysis reveals that, 
at 172 > 0, the interaction is marginally relevant. Usually 
the development of a strong coupling regime is accom- 
panied by a dynamical mass generation and the loss of 
conformal invariance at the strong coupling fixed point. 
This is indeed the case for the model of two marginally 
coupled S =1/2 spin chain with zigzag interaction 
We shall see, however, that due to chiral asymmetry of 
7i the effective interaction, as in the two-channel Kondo 
model, flows towards an intermediate fixed point where 
conformal invariance is recovered with a smaller central 
charge. 

The model (Q) is Bethe-Ansatz solvable |l6) (see also 
PJ). Here we present an exact solution for a U(l) ver- 
sion of the model, characterized by anisotropic interac- 
tion (52 — > 9\\i9-l)i which allows us to investigate the 
spectrum of the model and estimate asymptotics of the 
correlation functions. Our solution is based on a map- 
ping onto Majorana fermions. Using Abclian bosoniza- 
tion, we shall exploit the existence of a Toulouse-like 
point where the fermions are free We start with 

Abelian bosonization of the SU(2)i current Jo- Intro- 
ducing a massless bosonic field tp, we have (see, for in- 
stance, Appendix A of Ref. @): J^ R L = -|= d x ip R . L , 

J qrl = 2^T e Tl ^ VR - L . On the other hand, com- 
bining two Majorana fields, and £ 2 to form a Dirac 
fermion, (£ 2 + i^ 1 )/^/^, and then bosonizing it, we can 
express the SU(2)2 current I in terms of a bosonic field 
$ and the Majorana fermion £ 3 : Irl = d x <&R t L, 

ii T = ; 1 £% T KQ^ViTr® r,l An additional fermionic 

zero-mode operator n has been introduced to ensure the 
correct commutation relations. Then the Hamiltonian 7i 
can be written in the following bosonized form: 

n = v [{d xVR f + (d x p L f] + vi [{d x <s> R f + (d x <s> L ) 2 } 

- iy - e L d x e L ) 

+ -J- (dxPLdx^R + d x ip R d x $ L ) 
\/2tt 

2(vra ) 3/2 V > 

+ ^SrtCi^ cos (V4^$l + Vftetpii) ■ ( 5 ) 

2 (too) v ' 

We now perform a canonical transformation: 

if = cha $2 + sha $1, $ = cha $1 + sha $2 

§ — cha ©2 ^ sha 0i, = cha 0i — sha 02 (6) 
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where <d and 9 (respectively 81 and 82) are the dual 
fields associated with <p and $ (respectively $1 and $2). 
The cross terms d x ipd x Q in Eq. (^) can be eliminated by 
setting th2a = —g\\ /tt\/2 (vo + v\) . One immediately ob- 
serves that choosing tha = — -j= , which corresponds to a 
special (though nonuniversal) positive value of <?|j, (g 
4tt(vo + i?i)/3), the arguments of the two cosine terms 
in (||) become those of free fermions, cos (v / 47r$2L.fi) ■ 
Introducing a pair of Majorana fields, 77 and £, and us- 
ing the correspondence: ip RjL = (r] RtL + iCr,l)/V2 = 
(k/V2tto^) e ±i ^ F * 2 ' R ' i , we finally obtain: 



n = Y [(&$i) 2 + (^eo 2 ] - ^ {(Rd x ( R - (Ml] 

iu 2 



standard definitions of the currents of the three decou- 
pled chains. In this (ultraviolet) limit, one has a picture 
of three groups of independently propagating spinons. 
However, in the IR limit (\u>\ <C m), all Majorana bilin- 
ears in (H) are characterized by short-ranged correlations 
_ (since the Majorana fermions r\ and £ 3 are massive), im- 
plying that strongly fluctuating parts of the currents of 
individual chains are no longer independent; all of them 
contribute to the formation of a single, physical, current 
J z . In fact, the physical spinon represents a chirally 
asymmetric, strongly correlated state of three spinons. 
Consider, for instance, a right-moving -\/7r/2-kink of the 
field $1, representing a physical spinon with the spin pro- 
jection S z — 1/2. According to the exact relation 



-y [&d*& - ddtd] - -f [vRd x r, R - VLd xVL ] 



(7) 



Here m = g±/2irao, and the two renormalized velocities, 
Mi and U2, are expressed in terms of the surface and mid- 
dle chain velocities Vi and vq as u\ = (2v\ — vg)/3, u 2 = 
(2vq — «i)/3. The first two terms in Eq. (Q) describe 
completely decoupled free massless bosonic and Majo- 
rana fields, $1 and £, contributing to criticality with the 
central charge: C = 1 + 1/2 = 3/2. Therefore, at the new 
critical point Ti effectively represents a gapless S = 1/2 
spin chain and a critical Ising model. Coming back to 
the model (||) and adding the contribution of the singlet 
Majorana fermion £°, i.e. one more critical Ising model, 
the total central charge becomes Ct = 2. The remaining 
part of the Hamiltonian ((?]) has a spectral gap m and 
describes hybridization of the Majorana £ 3 and 77 fields 
with different chiralities. Since the canonical transforma- 
tion (^) does not mix Hi and Ti 2 , the Hamiltonian (0) 
still decouples into two commuting, chirally asymmetric 
parts. This reflects the chiral nature of the fixed point. 

Physical picture of elementary excitations. 
There are two different kinds of elementary excitations 
at the Toulouse point: magnetic excitation described by 
the field $1, and nonmagnetic, singlet, excitations asso- 
ciated with the two Majorana fermions £° and £. 

Notice that, due to the mixing of different degrees of 
freedom reflected in the canonical transformation (^) , the 
"physical" spinons, i.e. those defined as y/ir/2- kinks of 
the field $1 describing the effective 5 = 1/2 spin chain, 
should not be misleadingly identified as the spinons of 
the middle chain. To get a better understanding of the 
structure of spin excitations at the chiral fixed point, let 
us express the currents Jf , Jf and Jq in terms of the 
"physical" current J z = {1/ ^/2~K)d x $i. Using the trans- 
formation (||) at the Toulouse point, we find: 



J\,ir{l) = •Jr(L) _ g \P L (R)^ L (R) T £r(l)£r(l) ) J 

J 0R(L) = ~^L(R) + 1VR(L)(r(L)- (8) 

At energies \oj\ ^> m, where all Majorana fields can be 
considered as massless, Eqs. (0) transform back to the 



R(L) 



JlR(L) + J%R 



Jql(r) 



(9) 



following from (@), such an excitation is a combination 
of two right-moving spinons of the surface chains, each 
carrying the spin S z — 1/2, and a left-moving antispinon 
of the middle chain, with S z = —1/2. The rigidity of such 
a state is ensured by a finite mass gap in the {g—£, 3 ) sector 
of the model. This peculiar structure of the elementary 
spin excitations at the chiral fixed point is also reflected 
by the expression for the velocity u\. 

Apart from the nontrivial nature of the spinon, the 
chiral fixed point manifests itself in the existence of two 
new massless singlet excitations. The Majorana fermion 
£° describes collective excitations of singlet pairs formed 
on the two surface chains. The nature of the Majorana 
fermion £ is less transparent: it is a highly nonlocal object 
when expressed in terms of the original spin operators. 
The best way to understand the role of the singlet excita- 
tions is to combine the £ and C fields into a single Dirac 
fermion and then bosonize it. The corresponding mass- 
less bosonic field $ c resembles the scalar field describing 
the charge degrees of freedom in the Hubbard model away 
from half filling in the limit U = 00. We shall hence 
refer to the gapless Majorana fields as "pseudocharge" 
excitations which account for the central charge C = 1. 
The spin-pseudocharge separation is already manifest at 
the Toulouse point implying that the leading asymptotics 
of correlations functions will factorizc into the spin and 
pseudocharge contributions. 

Correlation functions. We shall now use the ex- 
act solution of the model, found at the Toulouse point, 
to calculate physical quantities of interest. From Eqs. 
(||) it follows that, as for a single S — 1/2 Heisen- 
berg spin chain, an external magnetic field H couples 
only to the massless field $1. The uniform suscep- 
tibility in units of g[iB is easily found to be x _1 = 
2ttMi. At low temperatures (T <C to) only the gap- 
less modes (£°, ( and $1) contribute to the specific 
heat. Using the general formula CV = ttCT/3v, we find: 
ttT (1 1 11 1 

6 \2 Vi 2 1*2 U\ / 
We have also computed all spin-spin correlation func- 
tions at the Toulouse point (details will be published else- 
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where |T3| ) . The leading asymptotics of the uniform part 
of the correlation functions coincide with that of the ef- 
fective S = 1/2 spin chain, in agreement with the above 
discussion. The nature of the chiral fixed point mani- 
fests itself in staggered correlations between the spins of 
the surface chains: 



(n a (x,t) n b (0,0)) 

Sab 



(x 2 + ufr 2 ) 1 / 2 (x 2 + vjr 2 ) 1 / 8 (x 2 + U 2 r 2)l/8 



(10) 



where a, b — (1,2). Notice the factorized contribution 
of the "pseudocharge" coming with the exponent 1/8. 
As seen from Eq. (|io|), there is an intrinsic velocity 
anisotropy in the pseudocharge sector which might be im- 
portant when considering the dynamical properties. We 
emphasize that the exponents in the correlation func- 
tions are universal and characterize a new universality 
class in spin ladders. This is the main result of our work. 
Another quantity of experimental interest is the NMR re- 
laxation rate 1 jT\ . It is not difficult to show that at low 
temperature, 1/Ti ~ VT, in contrast with the Heisen- 
berg chain where 1/Ti ~ const. 

Stability of the chiral fixed point. With all these 
results at hand, now we turn to the stability of the chiral 
fixed point. There are two important questions we shall 
adress. The first is related to the stability of the Toulouse 
point with gi = kept fixed, and the second is to examine 
the behavior of the system when one moves away from 
the point J± = 2 J x . 

We begin by stressing that our solution is stable pro- 
vided both u\ and ui are positive, the stability condi- 
tion thus being 1/2 < v§jv\ < 2. So the chiral fixed 
point is stable in a relatively broad range of velocities. 
When g» deviates from its Toulouse-point value, the 
Hamiltonian (]?]) picks up an extra term Sgi\ {d x $iL(RilR 
+d x ^iRCLilL)- Since the field r\ is massive, the expan- 
sion in 5g\\ does not introduce new infrared singularities, 
implying that the long-distance behavior of the correla- 
tion functions will not be modified except for a velocity 
and mass renormalization. Therefore, the solution at the 
Toulouse point captures all universal properties of the 
chiral fixed point and includes the case = g± where 
the Hamiltonian has the full SU(2) symmetry. 

Now we consider small deviations from the point g\ = 
0. At small enough \gi\ it is possible to investigate the 
effect of the backscattering term no • (ni + n2) as a weak 
perturbation to the chiral fixed point. This perturbation 
can be shown to be proportional to cos(v / 7r$ c ), which 
is a relevant operator with scaling dimension 1/4. We 
thus conclude that the backscattering operator opens a 
gap, A c , in the pseudocharge sector but has no effect on 
the magnetic (spinons) excitations. Standard scaling ar- 
guments give an estimate: A c ~ g 4 ^*' ' . The chiral fixed 
point is thus unstable in the far IR limit, and the sys- 
tem will flow to the C = 1 fixed point of the standard 
three-leg spin ladder. Of course, the very applicability 



of the perturbative approach to the chiral fixed point re- 
quires that A c -C m, the condition which can always be 
satisfied for sufficiently small g\. Under this condition, 
there exists an intermediate but still low-energy region 
A c <C E <C m where the C = 2 behavior caused by frus- 
tration is dominant. The physics in this region is univer- 
sal and cannot be understood without having recourse to 
the new chiral fixed point. At lower energies, E <C A c , 
the system will eventually cross over to the conventional 
critical C = 1 behavior. 

We think that this new fixed point discussed on the 
present paper might be responsible for new physics in 
many frustrated ladders such as the three-chain zigzag 
ladder. This model and the doped case are currently 
under study. We hope that the chiral-fluid critical state 
with all its physical properties will be observed in further 
experiments on spin-ladder systems. 
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